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1.  Introduction 


Currently,  a  common  practice  in  wind  modeling  for  small-arms  applications  is  to 
assume  that  a  constant  crosswind  exists  across  the  entire  range  to  a  target.  We 
know,  however,  that  wind  does  not  behave  this  way  in  reality.  Analysis  and  study 
of  the  behavior  of  wind  can  help  efforts  to  understand  and  reduce  the  delivery  error 
attributed  to  crosswind. 

The  first  goal  of  this  study  was  to  develop  a  methodology  to  generate  synthetic 
wind  fields  (i.e.,  the  crosswind  as  a  function  of  location  at  a  fixed  point  in  time) 
based  on  empirical  data.  The  ability  to  incorporate  synthetic  wind  fields  into 
existing  simulations  will  increase  those  models’  accuracy  and  provide  us  with  a 
better  understanding  of  weapon  accuracy.  The  differences  between  the  effects  of 
constant  crosswind  across  all  ranges  and  the  effects  of  variable  wind  are  likely  small 
at  close  ranges,  but  in  long-range  applications  the  difference  can  be  large,  and  we 
could  be  underestimating  or  overestimating  the  magnitude  of  deflection  due  to 
crosswind  at  range. 

This  could,  in  turn,  lead  to  an  overestimation  or  underestimation  of  the  gains  in 
accuracy  from  a  wind  estimation  and  correction  system.  Moreover,  the  ability  to 
model  variable  wind  will  give  us  the  ability  to  more  accurately  model  the 
perfonnance  of  these  systems  directly.  With  an  explicitly  modeled  synthetic  wind 
field  we  can  better  model  the  perfonnance  of  a  system  that  estimates  and  corrects 
for  that  wind. 

The  second  goal  of  this  study  was  to  develop  a  metric  called  Equivalent  Constant 
Crosswind  (ECC),  which  allows  for  the  effects  of  variable  crosswind  on  a 
downrange  projectile  to  be  modeled  without  having  to  explicitly  model  the 
variation  of  the  wind  at  every  step.  This  will  allow  for  easy  incorporation  of  wind 
fields  into  existing  models  without  extensive  changes  and  without  substantially 
increasing  the  runtimes  required. 

To  accomplish  the  study  goals,  we  analyzed  a  set  of  crosswind  data  obtained  from 
an  anemometer  array  at  Fort  Benning  in  December  2012  and  investigated  several 
methodologies  to  simulate  the  observed  variation  of  crosswind  with  range, 
including  treating  each  anemometer  independently,  autoregressive  models  of  order 
p  (where  the  previous  p  crosswind  measurements  are  used  to  detennine  the  next 
crosswind  value),  and  sampling  against  a  multivariate  normal  distribution. 

We  ran  a  Monte  Carlo  analysis  over  a  matrix  of  input  parameters  to  determine  the 
effect  of  those  parameters  on  ECC.  The  parameters  most  relevant  to  this  analysis 
were  shooter-to-target  range,  muzzle  velocity,  projectile  mass,  drag  coefficient 
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exponent,  wind  standard  deviation  (i.e.,  the  standard  deviation  of  the  crosswind 
magnitude  at  the  measurement  points  over  time),  and  spatial  wind  correlation. 

Section  2  will  summarize  the  analysis  done  on  the  data  set  obtained  at  Fort  Benning 
in  2012.  Section  3  will  detail  the  methodology  for  using  McCoy’s  lag  time 
methodology  alongside  the  Weinacht-Cooper-Newill  (WCN)  projectile  flyout 
methodology  to  obtain  measures  of  ECC.  Section  4  will  briefly  summarize  the  ECC 
results  obtained  from  the  observed  data  and  the  assumptions  made  in  obtaining 
them.  Section  5  will  outline  the  different  approaches  used  to  model  the  observed 
wind  behavior.  Section  6  will  present  the  recommendations  for  how  to  use  the 
model  that  fits  the  data  best  along  with  the  effects  of  changing  different  variables. 
Section  7  will  present  the  conclusion. 

2.  Analyzing  Wind  Data 

In  December  2012,  a  set  of  wind  values  was  collected  at  Fort  Benning.  Ten 
anemometers  (Ai-Aio)  were  placed  along  a  straight  line  spaced  110  m  apart.  Wind 
speed  and  direction  were  recorded  every  second  for  just  over  3  h,  and  from  these 
values  the  crosswind  and  range  wind  were  derived.  For  our  purposes,  crosswind  is 
defined  as  the  component  of  the  wind  perpendicular  to  the  line  connecting  the 
anemometers,  with  positive  crosswind  values  to  the  right  as  you  look  from  Ai  to 
A  io,  and  range  wind  is  defined  as  the  component  of  the  wind  parallel  to  the  line 
connecting  the  anemometers,  with  positive  range  wind  in  the  direction  of  ascending 
anemometer  numbers.  Because  crosswind  is  a  much  bigger  contributor  to 
deflection  at  range  than  range  wind,  we  only  analyzed  the  crosswind  component  of 
the  data  set. 

The  crosswind  measurements  from  each  anemometer  were  analyzed  for  normality 
and  were  found  in  most  cases  to  be  close  to  normal.1  A  couple  of  the  locations 
showed  a  somewhat  pronounced  skew.  Figure  1  is  a  histogram  of  the  measured 
crosswind  velocities  at  A3  that  appear  to  be  roughly  normal.  Figure  2  is  a  histogram 
of  the  measured  crosswind  velocities  at  A4  that  show  a  pronounced  skew.  To 
simplify  calculations  when  modeling  variable  wind  fields,  the  wind  measurements 
at  each  station  were  assumed  to  be  normal.2 
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Fig.  1  Histogram  of  measured  crosswind  values  from  A3 
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Fig.  2  Histogram  of  measured  crosswind  values  from  A4 

Mean  and  standard  deviations  of  the  crosswind  velocities  at  each  anemometer  were 
also  obtained  from  the  data  and  are  shown  in  Table  1  and  graphically  in  Figs.  3  and 
4.  For  these  calculations,  each  measurement  (taken  at  the  same  location 
approximately  once  per  second)  was  assumed  to  be  an  independent  sample. 
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Table  1  Crosswind  means  and  standard  deviations  (in  meters/second) 


Ai 

A2 

As 

A4 

As 

A6 

A7 

A8 

A9 

A10 

Mean 

-1.55 

-1.59 

-1.71 

-1.77 

-1.47 

-2.69 

-1.41 

-2.74 

-2.93 

-3.01 

Std.  Dev. 

1.28 

1.55 

1.34 

1.59 

1.58 

1.99 

1.61 

1.83 

1.59 

1.72 

Mean  Crosswind  Velocity 


Fig.  3  Mean  crosswind  velocity  vs.  range 


Crosswind  Velocity  Standard  Deviation 


Fig.  4  Crosswind  standard  deviation  vs.  range 

For  the  remainder  of  the  report,  the  mean  and  standard  deviation  of  the  crosswind 
at  anemometer  An  will  be  denoted  pn  and  an,  respectively.  Instances  of  the 
crosswind  at  anemometer  An  will  be  denoted  Wn. 

To  model  the  behavior  of  the  crosswind  accurately  we  first  had  to  determine  the 
correlation  between  the  observed  wind  values.  We  can  easily  determine  the 
correlation  coefficient  p  between  stations  An  and  Am  by  the  following  equation 
(Box  1976): 
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(1) 


P(Ani  Am) 


cov(Wn,Wm ) 
anam 


This  generated  Table  2,  which  lists  correlation  coefficients  between  each  pair  of 
stations.  A  correlation  coefficient  of  1 .0  indicates  total  positive  correlation  between 
2  variables,  meaning  that  a  linear  equation  describes  their  relationship  perfectly.  A 
correlation  coefficient  of  0  indicates  no  linear  correlation  between  2  variables. 
Values  in  between  indicate  varying  levels  of  positive  correlation  between  variables. 


Table  2  Pairwise  correlation  coefficients  between  anemometers 


Ai 

A2 

As 

A4 

As 

A6 

A7 

A8 

A9 

A10 

Ai 

1.00 

0.27 

0.16 

0.21 

0.20 

0.25 

0.21 

0.17 

0.17 

0.18 

^2 

0.27 

1.00 

0.19 

0.20 

0.29 

0.23 

0.22 

0.19 

0.18 

0.17 

^3 

0.16 

0.19 

1.00 

0.23 

0.27 

0.30 

0.24 

0.21 

0.21 

0.21 

a4 

0.21 

0.20 

0.23 

1.00 

0.32 

0.32 

0.31 

0.28 

0.25 

0.25 

a5 

0.20 

0.29 

0.27 

0.32 

1.00 

0.36 

0.34 

0.30 

0.29 

0.24 

a6 

0.25 

0.23 

0.30 

0.32 

0.36 

1.00 

0.53 

0.43 

0.39 

0.35 

Ay 

0.21 

0.22 

0.24 

0.31 

0.34 

0.53 

1.00 

0.48 

0.36 

0.34 

As 

0.17 

0.19 

0.21 

0.28 

0.30 

0.43 

0.48 

1.00 

0.43 

0.41 

a9 

0.17 

0.18 

0.21 

0.25 

0.29 

0.39 

0.36 

0.43 

1.00 

0.47 

A10 

0.18 

0.17 

0.21 

0.25 

0.24 

0.35 

0.34 

0.41 

0.47 

1.00 

For  the  remainder  of  this  report,  correlation  coefficients  between  stations  An  and 
Am  will  be  denoted  pn  rn.  There  is  no  difference  between  pn  m  and  pmn,  but  the 
preference  will  be  for  the  smaller  number  to  be  listed  in  the  subscript  first. 

In  the  next  section,  we  will  explain  in  some  depth  the  methodology  used  to  analyze 
the  observed  crosswind  data  and  ultimately  to  generate  synthetic  wind  fields. 


3.  Calculation  of  Horizontal  Deflection  and  Equivalent  Constant 
Crosswind 

To  determine  the  horizontal  deflection  imparted  by  a  variable  crosswind  on  a 
projectile  at  its  target,  we  used  a  combination  of  equations  detailing  deflection  due 
to  variable  crosswind  (McCoy  1976)  and  the  WCN  equations  detailing  trajectories 
of  direct- fire  munitions  (Weinacht  et  al.  2005). 

3.1  McCoy  Methodology  for  Deflection  at  Range  Due  to 
Crosswind 

To  determine  the  horizontal  deflection  (ZR.)  imparted  at  range  R  by  a  constant 
crosswind  occurring  between  ranges  X^x  andAj  (for  Xi,Xi+1  <  R),  McCoy  (1976) 
gives  the  following  equation3: 
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where 


J Rt 


W; 


(i— i,i)  [ 


tR  tXi_±  — 


R-Xi _ 


Vx, 


*R  tx. 


R-Xi 


VXi 


(2) 


WZ(i_ i0  =  magnitude  of  crosswind  between  ranges  Xt_t  and  Xt, 
tx i  =  time  of  flight  of  the  projectile  to  range  Xt,  and 
Vx.  =  residual  velocity  of  the  projectile  at  range  XL. 

R _ X  ■ 

The  tenn  tR  —  tx.  — -^—L  in  Eq.  2  represents  the  difference  between  actual  flight 

time  between  range  Xt  and  range  R  and  the  time  between  the  same  ranges  in  a 
vacuum.  This  is  commonly  referred  to  as  “lag  time”. 

If  there  is  a  set  of  monotonically  increasing  ranges  {X0,X1,X2, ...  ,Xn],  < 

R  V  i  <  n,  and  a  corresponding  set  of  crosswinds  {^z(04).  wzil}2y  -  -  wz^n_ln)} 
where  crosswind  between  ranges  X[_1  and  Xi  is  constant  and  equal  to  Wz  _t  ,  then 
the  total  horizontal  deflection  at  range  R  due  to  crosswind  is  equal  to 


It  stands  to  reason  that  there  exists  a  constant  crosswind  that  would  have  produced 
the  same  horizontal  deflection  as  any  particular  variable  wind  field.  If  we  apply 
Eq.  2  across  a  trajectory  with  constant  crosswind  Wz,  range  R,  and  time  of  flight  t, 
we  arrive  at  the  following  equation  for  ZR : 


ZR  =  Wz 


(4) 


By  combining  Eqs.  3  and  4  and  solving  for  Wz,  the  constant  crosswind  that  would 
produce  the  same  deflection,  or  ECC,  is  given  by 


ECC  =  Wz 


(5) 


3.2  Incorporation  with  WCN  Model 

The  calculation  of  ECC  depends  on  time  of  flight  and  residual  velocity  for  every 
range  Xt  used  as  an  input.  To  determine  these  quantities  at  a  given  range  for  a 
general  case  that  correctly  models  flight  in  the  transonic  and  subsonic  regimes,  we 
use  the  WCN  methodology  (Weinacht  2014).  A  closed-form  solution  can  be  used 
if  the  projectile’s  velocity  remains  supersonic.  This  solution  is  outlined  in 
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Appendix  A.  For  cases  where  flight  crosses  velocity  regimes,  we  use  Oberle’s 
implementation  of  the  WCN  methodology  (Oberle  2009). 

Under  the  assumption  that  the  drag  curve  is  continuous,  Oberle’s  model  requires 
the  following  inputs: 

•  Muzzle  velocity 

•  Drag  coefficient  at  the  muzzle  velocity 

•  Drag  coefficient  exponent  in  the  supersonic  velocity  regime 

•  Transition  Mach  number  between  transonic  and  supersonic  velocity 
regimes 

•  Transition  Mach  number  between  subsonic  and  transonic  velocity  regimes 

•  Drag  coefficient  in  the  subsonic  velocity  regime  (CD  \Suh) 

•  Muzzle  retardation 

If  we  are  to  investigate  notional  projectiles  with  modified  mass  and  velocity,  we 
must  have  expressions  that  allow  us  to  find  the  corresponding  muzzle  retardation 
and  drag  coefficient  at  the  muzzle  velocity.  Muzzle  retardation  is  expressed  as 
follows: 

(%)0  =  ^pV°s"'c°'v°-  <6) 

where 


m  =  projectile  mass, 
p  =  air  density, 

Sref  =  reference  area  of  the  projectile,  and 
CD\v0  =  drag  coefficient  at  muzzle  velocity. 


If  we  know  the  air  density  as  well  as  the  mass,  muzzle  retardation,  and  reference 
area  of  the  projectile,  we  can  solve  Eq.  6  for  Cd\Vq  as  a  function  of  known  quantities 

and  not  specify  it  explicitly. 


CD  | 


Vo 


pVoSref 


(7) 


Oberle  includes  the  following  expression  for  deriving  muzzle  retardation  from  a 
known  set  of  values  for  bullets  with  the  same  basic  shape  and  reference  area  but 
with  different  mass  and  muzzle  velocity  values. 
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(8) 


/‘f  l7)  _  f dv  )  f  m Known  \  f  V New 

\ds  J New  \dsj Nnown  '  mNew  '  W Known 

Quantities  with  the  Known  subscript  refer  to  known  values  and  quantities  with  the 
New  subscript  refer  to  modified  velocity  and  mass.  This  means  that  if  there  is  an 
array  of  associated  muzzle  velocity,  muzzle  retardation,  and  mass  and  drag 
coefficient  exponent  values,  then  muzzle  retardation  can  be  expressed  as  a  function 
of  those  values  and  modified  mass  and/or  muzzle  velocity  values. 

Thus  muzzle  retardation  and  drag  coefficient  at  the  muzzle  can  be  removed  as 
independent  variables  in  favor  of  mass,  air  density,  and  reference  area,  and  the  list 
of  required  inputs  is  modified  to  the  following: 

•  Muzzle  velocity 

•  Projectile  mass 

•  Projectile  reference  area/diameter 

•  Air  density 

•  Set  of  associated  supersonic  drag  coefficient  exponents,  muzzle  velocities, 
muzzle  retardations,  and  masses 

•  Transition  Mach  number  between  transonic  and  supersonic  velocity 
regimes  (assumed  to  be  1 .05  for  this  analysis) 

•  Transition  Mach  number  between  subsonic  and  transonic  velocity  regimes 
(assumed  to  be  0.85  for  this  analysis) 

•  Drag  coefficient  in  the  subsonic  velocity  regime 

Using  these  inputs,  the  trajectories  can  then  be  run  using  Oberle’s  implementation 
of  the  WCN  model. 

4.  Analysis  of  Equivalent  Constant  Crosswind  from  Wind  Data 

We  first  analyzed  the  behavior  of  the  ECC  of  the  Fort  Benning  wind  data.  The 
results  of  this  analysis  became  the  standard  against  which  we  judged  the 
performance  of  our  proposed  synthetic  wind  field  models.  The  data  set  included 
crosswind  values  at  each  of  the  10  anemometers  at  1 1,085  different  fixed  times.  We 
treated  these  as  independent  samples  for  the  purpose  of  this  analysis.  We  assumed 
that  the  crosswind  values  measured  at  each  station  at  the  time  the  shot  was  fired 
would  be  the  crosswind  values  experienced  by  the  projectile  when  it  arrived  at  that 
station.  That  is,  time  variation  of  the  crosswind  during  a  trajectory  did  not  enter  into 
the  calculation  of  deflection.  Second,  we  assumed  that  the  crosswind  at  any 
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trajectory  location  was  equal  to  the  measured  crosswind  value  at  the  closest 
anemometer. 


Adjacent  anemometers  were  separated  by  1 10  m  in  a  straight  line  from  0  to  990  m. 
We  determined  that  our  ranges  of  interest  were  every  55  m  from  the  shooter  to  the 
target;  that  is,  at  each  anemometer  and  at  each  point  equidistant  from  2  adjacent 
anemometers.  This  split  the  range  into  18  subsections  with  each  assumed  to  have 
constant  wind  within  them.  To  find  the  ECC  for  a  given  range  R,  we  used  Eq.  2  to 
determine  the  deflection  induced  at  R  by  the  wind  on  each  interval  and  then  we 
summed  those  deflections  and  input  the  sum  into  Eq.  5  to  obtain  the  ECC. 


We  used  the  WCN  methodology  to  find  the  time  of  flight  and  residual  velocity  at 
each  range  of  interest  using  the  M855/M4  weapon/ammunition  combination 
V0  —  866  m/s,  m  —  4  g,  d  —  5.56 mm,  (  — j  =  —1.028,  n  —  0.53,  CD\Sub  — 

\ds/Q 

0.198,  p  =  1.225  kg/mA3.  ECC  depends  on  the  variables  that  describe  the 
behavior  of  the  wind  (mean  and  standard  deviation  at  each  location  and  correlation 
between  locations)  as  well  as  the  inputs  that  describe  the  flight  of  the  particular 
projectile  being  analyzed.  Section  6  and  Appendix  B  go  into  greater  detail  about 
how  much  the  ECC’s  behavior  depends  on  each  of  these  inputs.  Figure  5  shows  the 
mean  ECC  value  versus  range  from  the  Fort  Benning  data  set  (/r EccB )>  and  Fig.  6 
shows  the  standard  deviation  of  ECC  versus  range  (<7ECCb). 


Range  (m) 


Fig.  5  Mean  values  of  HECCg  vs.  range  observed  at  Fort  Benning 
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Fig.  6  Standard  deviation  values  of  < tECCb  vs.  range  observed  at  Fort  Benning 

The  value  of  \iECcb  decreases  with  range  because  the  means  of  the  winds  get 
stronger  in  the  negative  direction  as  range  increases,  as  shown  in  Table  1 .  The  value 
of  oECcb  decreases  with  range  until  halfway  downrange,  at  which  point  it  begins  to 
increase  gently.  The  increase  in  the  longer  ranges  appears  to  be  a  consequence  of 
standard  deviation  of  measured  crosswind  increasing  with  range,  shown  in  Table  1, 
and  the  increase  in  correlation  among  anemometers  A6-A10,  shown  in  Table  2  (the 
correlation  values  among  anemometers  A1-A5  range  from  0.16  to  0.32;  the 
correlation  values  among  anemometers  A6-A10  range  from  0.34  to  0.53). 

The  distributions  of  ECC  at  each  range  were  found  to  behave  similarly  to  the 
distributions  of  the  crosswind  measurements  at  each  range.4  Figure  7  shows  the 
distribution  of  ECC  values  at  110  m,  which  appears  to  be  approximately  normal. 
Figure  8  shows  the  distribution  of  ECC  values  at  990  m,  which  has  some 
irregularity  including  a  small  secondary  local  maximum  to  the  left  of  the  mean. 
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Fig.  7  Histogram  of  calculated  ECC  values  at  110  m 


ECC  at  990m 

700 
600 
500 

> 

|  400 
S'  300 

s_ 

LI- 

200 

100  I 

q  _ ^rnfflT  I 

LDLncnmi^THLncnror^THLncTiror^THLnCTir^rocnLnTHr^ro 

LnoOTHTHrgrgrgroro«^-^f^LnLDtDii3ii3rgrgT-iTHTHOO 

^t-org-sj-iDoooiN'sj-iDoooiN^-tDooorsiLnmvHair^Lnm 

^r^oroiuairoiucniNLncniNLnoovHLnoooOLniNoOLnrgai 

‘I'oi^oi^jiLnqLnriiOriivNivmcismqLnq^oi'JM 

KKupvDLOLnLn'^^-rofrirlrlrHrHOOOOT-iT-iT-irlrl 

ECC  (m/s) 

Fig.  8  Histogram  of  calculated  ECC  values  at  990  m 

5.  Generating  Synthetic  Wind  Fields 

Next  we  attempted  to  find  a  model  that  could  reproduce  the  behavior  of  the  wind 
in  the  Fort  Benning  data  set  as  closely  as  possible.  We  first  generated  a  wind  field 
by  drawing  the  wind  at  each  station  independently  using  /i„  and  on  at  each  station 
An.  This  produced  oECC  results  that  diverged  unacceptably  from  (JECCb  ■  Next  we 
implemented  an  autoregressive  model,  which  generated  each  wind  measurement 
from  some  number  of  the  measurements  that  preceded  it.  This  model  generated 
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°ecc  results  much  closer  to  ctECCb  ,  but  the  fact  that  it  depended  on  assigning  an 
arbitrary  directionality  to  the  wind  measurement  locations  introduced  an  ambiguity 
to  its  implementation  that  led  us  to  investigate  alternate  methodologies.  Further 
discussion  of  the  foregoing  methodologies  is  included  in  Appendix  C.  To  refine  the 
results  from  the  independent  and  autoregressive  models,  it  is  necessary  to  find  a 
model  that  not  only  preserves  the  distributions  of  measured  wind  at  each  station, 
but  the  correlations  between  all  pairs  of  stations. 

Upon  investigation,  a  multivariate  normal  distribution  addresses  all  of  these 
concerns.  If  we  let  1)  W  be  an  n-dimensional  vector  whose  components  are  random 

variables  W1,  W2, ... ,  Wn,  2)  X,  the  covariance  matrix,  be  an  n  x  n  matrix  whose  i,jth 

— ) 

entry  is  the  covariance  of  WL  and  Wj,  3)  Z  be  an  n-dimensional  vector  whose 
components  are  independent  standard  normal  variates,  and  4)  p  be  an 
n-dimensional  vector  whose  ith  entry  is  E  (I/F,),  then 

W  =  p  +  AZ  ,  (9) 

where  A  is  an  n  x  n  matrix  such  that  A  A1'  —  X  (Gentle  2009). 

This  is  a  comparatively  easy  model  to  implement.  The  most  complicated  part  is  the 
calculation  of  the  matrix  A.  One  type  of  matrix  with  this  property  is  known  as  the 
Cholesky  decomposition,  and  there  is  a  straightforward  process  by  which  it  can  be 
calculated  (Golub  1996). 


The  number  of  iterations  to  run  in  each  Monte  Carlo  simulation  was  determined  by 
the  following  fonnula  (Oberle  2015): 


n  = 


(10) 


where 


n  =  number  of  iterations, 

zay  =  standard  normal  z-score  corresponding  to  a  confidence  interval  of 

100(1- a)%, 

o  =  predicted  standard  deviation  of  the  data,  and 
A  =  half  length  of  the  confidence  interval. 

We  calculated  the  number  of  iterations  using  the  chosen  value  of  wind  standard 
deviation  as  er,  a  confidence  level  100(1-  a)%  of  98%,  and  an  interval  half-length 
A  of  0.01  m/s.  We  used  a  wind  standard  deviation  of  1.30  m/s  to  yield  91,461 
iterations,  which  we  rounded  up  to  100,000. 
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With  100,000  runs,  it  is  easy  to  verify  that  nn  and  ern  is  preserved  for  all  n  and  that 
all  of  the  pairwise  correlations  between  stations  are  preserved.  Figure  9  shows  the 
results  of  a  comparison  of  (tECCmvn  with  oECCb. 
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Fig.  9  Multivariate  normal  sampling  oECC  results 

Using  this  method,  oECCmvn  differed  from  oECCb  by  no  more  than  1%  at  any  range. 
As  with  the  other  models,  [J-eccmvn  remains  within  1%  of  [J-eccb  at  all  ranges. 

6.  Effects  of  Wind  Field  Changes 

The  aim  of  this  report  is  to  develop  a  model  that  can  be  used  to  generate  general 
crosswind  fields  based  on  a  set  of  inputs.  We  can  make  small  changes  to  the 
observed  data  to  see  how  changes  to  various  inputs  affect  the  ECC  output. 

Changing  the  wind  standard  deviation  roughly  represents  a  vertical  stretch  or 
compression  of  the  aECC  versus  range  curve,  as  shown  in  Fig.  10.  Uniformly  halving 
an  for  all  n  has  the  effect  of  halving  oECc  at  each  range.  Adding  50%  to  on  for  all 
n  has  the  effect  of  adding  50%  to  aECC  at  each  range. 
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Fig.  10  Comparison  of  oECC  at  different  values  of  wind  standard  deviation 


Changing  the  correlation  values  changes  the  slope  of  the  oECC  versus  range  curve. 
In  Fig.  1 1  we  tested  3  levels  of  correlation.  The  baseline  case  was  the  one  outlined 
in  Table  2.  In  the  low  correlation  case  we  subtracted  0.1  from  each  observed 
correlation.  In  the  high  correlation  case  we  added  0.1  to  each  observed  correlation 
(apart  from  the  self-correlations,  which  were  left  at  1 .0).  The  higher  the  correlation, 
the  closer  oECC  stays  to  the  initial  chosen  value  of  the  wind  standard  deviation. 


Range  (m) 


■  High  Correlation 


■  Regular  Correlation 


Low  Correlation 


Fig.  11  Comparison  of  oECC  at  different  levels  of  correlation 


Approved  for  public  release;  distribution  is  unlimited. 

14 


In  the  data  set  obtained  from  Fort  Benning,  oy  <  on  for  all  n  A  1,  yet  aECC  <  <j1 
for  all  n  >  1,  despite  a  small  increase  starting  around  500  m.  This  is  because  the 
winds  at  individual  anemometers  are  only  weakly  correlated,  meaning  there  is  a 
strong  possibility  that  winds  at  separate  stations  act  in  opposing  directions  and 
partially  cancel  each  other  out.  As  shown  in  Fig.  1 1 ,  with  higher  correlations  among 
anemometers,  oECC  falls  less  with  range  than  it  does  with  lower  correlations  since 
more  highly  correlated  winds  are  more  likely  to  act  in  the  same  direction. 

The  small  increase  in  oECC  observed  from  500  to  1,000  m  with  this  particular  data 
set  is  caused  by  2  factors:  1)  on  tends  to  increase  with  range  and  2)  pnm  is 
uniformly  greater  when  6  <  n,  m  <  10  than  when  1  <  n,  m  <  5.  Figure  1 1  shows 
that  an  increase  in  correlation  causes  an  increase  in  <rECC.  In  this  case,  when  a 
projectile  flies  through  the  ranges  covered  by  anemometers  A6-A10,  it  is  more  likely 
to  encounter  crosswinds  acting  in  concert  with  each  other  and  therefore 
compounding  the  deflection  at  the  target  rather  than  canceling  each  other  out. 

We  found  that  aECC  depended  on  an  and  pnm  but  not  pn.  Conversely,  pECC 
depended  only  on  pn  and  not  on  on  or  pn7n.  Figure  12  shows  \iECC  for  the  3 
correlation  scenarios  previously  described,  and  the  curves  overlap  completely.  The 
result  is  the  same  for  any  cases  that  use  the  same  pn  for  all  n. 
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Fig.  12  Comparison  of  nECC  at  different  levels  of  correlation 

Changes  in  ballistic  inputs  had  very  little  impact  on  the  value  of  oECC •  However, 
this  does  not  mean  the  deflection  produced  by  crosswind  is  almost  identical  for  all 
rounds  irrespective  of  drag,  muzzle  velocity,  and  mass.  Rounds  with  higher  muzzle 
velocities  as  well  as  heavier  rounds  experience  less  deflection  due  to  crosswind 
because  their  time  of  flight  to  the  target  is  lower.  The  drag  exponent  has  a  second- 
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order  effect  on  the  round’s  velocity  over  its  flight  and  thus  plays  a  much  smaller 
role  than  mass  and  velocity  in  determining  the  deflection  of  the  round  due  to 
crosswind.  Graphs  presenting  these  results  are  included  in  Appendix  C.  The  fact 
that  the  ballistic  inputs  had  little  impact  on  the  aECC  values  indicates  that  we  can 
get  a  good  approximation  of  aECC  for  many  classes  of  projectiles  by  explicitly 
modeling  it  for  only  one  configuration. 

7.  Conclusion 


This  report  has  presented  a  way  to  generate  synthetic  wind  fields  over  ranges 
relevant  to  small-arms  fire  and  argues  that  variable  wind  fields  appear  to  behave 
according  to  a  multivariate  normal  distribution.  First  we  tried  modeling  wind 
measurements  independently,  but  the  modeled  wind  behavior  did  not  align  well 
enough  with  the  observed  data.  An  autoregressive  model  brought  us  closer  to 
matching  the  observed  results  but  exhibited  ambiguity  about  exactly  how  to 
configure  the  model.  In  addition,  the  autoregressive  model  used  only  wind 
measurements  from  one  direction  to  generate  new  wind  measurements,  but  ideally 
the  model  should  be  independent  of  the  direction  the  shooter  is  facing.  Sampling 
against  a  multivariate  normal  distribution  was  found  to  strongly  correlate  to  the 
observed  wind  across  all  metrics  and  operates  independently  of  the  locations  of  the 
shooter  and  target. 

We  developed  a  metric  to  measure  the  effect  induced  by  a  variable  crosswind 
experienced  by  a  projectile  throughout  its  trajectory.  For  any  deflection  induced  by 
crosswind  at  the  target  range  there  has  to  be  a  constant  crosswind  that,  if  applied 
over  the  whole  range,  would  have  produced  that  same  deflection.  This  was 
designated  the  ECC,  and  analysis  was  performed  to  characterize  its  behavior  with 
respect  to  a  number  of  different  inputs  to  the  model. 

We  found  that  the  ECC  at  a  given  range  behaved  according  to  a  normal  distribution, 
and  that  <rECC  tends  to  fall  as  range  increases  for  cases  where  the  wind  standard 
deviation  and  correlation  are  constant  across  the  whole  range.  In  the  case  of  the 
data  from  Fort  Benning,  the  wind  standard  deviation  varied  at  each  station,  which 
along  with  a  variation  in  correlation  led  to  a  small  increase  in  oECC  at  longer  ranges. 
Even  in  the  case  of  the  data  from  Fort  Benning,  though  the  value  of  oECC  increases, 
the  maximum  value  it  attains  is  at  least  20%  less  than  the  minimum  wind  standard 
deviation  measured  at  any  single  location.  This  implies  that  using  a  constant 
standard  deviation  derived  from  any  one  measurement  location  to  model  crosswind 
across  the  whole  range  will  tend  to  overestimate  the  horizontal  error  due  to 
crosswind,  though  the  overestimation  is  comparatively  small  at  shorter  ranges.  The 
factors  with  the  greatest  impact  on  aECC  were  the  range  of  projectile  travel,  the 
standard  deviation  of  crosswind  measurement,  and  correlation  between  crosswind 
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measurements.  Ballistic  characteristics  of  the  projectile  such  as  mass,  drag 
coefficient  exponent,  and  muzzle  velocity  had  comparatively  little  impact  on  the 
behavior  of  ECC. 

The  mechanism  driving  the  decline  in  oECC  with  range  appears  to  be  decreasing 
correlation  of  crosswind  measurements  with  increasing  spatial  separation.  Any  2 
winds  separated  by  only  a  short  distance  tend  to  be  more  strongly  correlated,  but  as 
range  between  crosswind  measurements  increases,  that  correlation  fades.  The 
crosswind  values  at  distant  ranges  become  only  weakly  correlated  with  the 
crosswind  value  at  the  shooter’s  location,  which  makes  downrange  winds  more 
likely  to  push  in  opposing  directions  the  farther  the  projectile  flies. 

The  ability  to  model  synthetic  wind  fields  will  better  allow  analysis  of  the 
effectiveness  of  wind  correction  methodologies  where  wind  is  read  at  multiple 
points  downrange  compared  with  the  effectiveness  of  correcting  using  the  wind  at 
the  shooter’s  location  in  the  presence  of  variable  wind.  Further  experiments  could 
provide  better  data  for  long-tenn  wind  standard  deviation  and  spatial  wind 
variability  at  various  locations,  times  of  day,  and  times  of  year. 
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8.  Notes 


1 .  Anderson-Darling  and  Kolmogorov-Smimov  tests  indicated  a  departure  from 
nonnality  but  show  very  high  sensitivity  for  sample  sizes  as  large  as  those  used 
here;  skewness  test  indicated  approximate  symmetry  with  a  persistent  small- 
to-moderate  negative  skew  at  all  ranges,  and  kurtosis  test  indicated  all  ranges 
were  slightly  leptokurtic. 

2.  Time  variation  of  wind  was  not  taken  into  account  for  this  analysis;  further 
study  of  variation  of  wind  in  both  space  and  time  may  refine  understanding  of 
the  underlying  distribution  of  the  wind. 

3.  Notation  has  been  adjusted  from  what  appears  in  the  source. 

4.  As  before,  Anderson-Darling  and  Kolmogorov-Smirnov  indicated  departures 
from  normality;  skewness  indicated  approximate  symmetry  with  a  small  but 
persistent  negative  skew  across  all  ranges  and  kurtosis  indicated  that  all  but  the 
2  furthest  ranges  were  somewhat  leptokurtic. 
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Appendix  A.  Synthesis  of  McCoy  and  Weinacht-Cooper-Newill 

Methodologies 
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Weinacht1  gives  the  following  expressions  to  detennine  the  time  of  flight  and 
residual  velocity  of  a  projectile  in  the  supersonic  range: 


1 


.  j  ,j  .  fdV\  XAn 
Vy.  =  Vn  11  +  n  —  —\ 

Xl  \dsJ0V0) 

,  n  A  0  . 

(A-l) 

(dV\ 

Vx.  =  V0e^dsJov0 

,  n  —  0  . 

(A-2) 

1  f„  ,  (dv\  XA1 

txt  =  f  ^ fdv A  1  +n(dJ  v 

Ll  Kds'°Vo) 

1 

n  —  1  ,  n  =£  0,n  =£  1 . 

(A-3) 

1  .  (.  ,  fdV\  XL' 

txt  =  (dv^  lni1  +  LJ  „ 

}  ,n  =  1  • 

(A-4) 

!  r  A] 

txi  =  (dV]  1  e  UJ°y° 
Ids  Jo  L  J 

,  =  0 , 

(A- 5) 

where 

E0  =  muzzle  velocity, 

( — )  =  muzzle  retardation  (s'1),  and 

\ds/ 0 

n  =  exponent  defining  the  shape  of  the  drag  versus  Mach  number  curve  (unitless). 

These  equations  can  be  combined  with  Eq.  2,  given  by  McCoy,2  to  produce  closed- 
form  expressions  for  the  horizontal  deflection  imparted  at  range  R  on  a  projectile 
that  remains  supersonic  throughout  its  flight  by  a  constant  crosswind  between 
ranges  X[_1  and  Xi  in  terms  of  crosswind  velocity  between  j  and  Xh  muzzle 
velocity,  muzzle  retardation,  and  the  exponent  defining  the  shape  of  the  drag  curve. 


Z„.  = 


wz(i-u)  r°+(f)0((n-i)fl+^ 


zRi  =  wz 


^-1)(a  [  (1+n(£)ogf  J 

w  {  1  in  Vo+^>oXi  | 

^  l(f  )„  M§ )„*- J  Mf )/.  Mf )/t-r 


,  n  ^  0,  n  ^  1  .  (A-6) 


n  =  1 .  (A-7) 


Zp.  =  WZr  , 


1-°  (Kf)0 


M  +  (*s\vlo  -  e  ids) 0vo 


,n  =  0  .  (A-8) 


'Weinacht  PA.  Direct-fire  trajectory  model  for  supersonic,  transonic,  and  subsonic  projectile  flight. 
Aberdeen  Proving  Ground  (MD):  Army  Research  Laboratory  (US);  2014  July.  Report  No.:  ARL-TR-6998. 
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Appendix  B.  Supplementary  Monte  Carlo  Results 
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Figures  B-l  through  B-4  show  the  standard  deviation  of  the  Equivalent  Constant 
Crosswind  (<JEcc)  versus  range  and  examining  the  changes  to  oECC  with  changes  in 
drag  coefficient  exponent,  muzzle  velocity,  projectile  mass,  and  subsonic  drag 
coefficient.  The  drag  coefficient  exponents  and  masses  were  chosen  from  among  a 
list  of  fielded  small-arms  projectiles  and  thus  are  not  evenly  spaced. 

Changing  the  drag  coefficient  exponent  alone  has  very  little  effect  on  oECC 
(approximately  1%)  at  any  range,  as  shown  in  Fig.  B-l. 
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Fig.  B-l  Comparison  of  different  values  of  n 

Changing  the  projectile’s  muzzle  velocity  has  more  of  an  effect  on  the  value  of 
<jECc  than  changing  the  drag  coefficient  exponent  (up  to  5%  in  the  most  extreme 
cases  examined),  but  that  effect  is  seen  mostly  at  extended  ranges  (Fig.  B-2).  Even 
at  the  largest,  though,  the  effects  on  aECC  induced  by  increasing  or  decreasing 
muzzle  velocity  are  dwarfed  by  the  effects  of  changing  wind  standard  deviation  and 
spatial  correlation. 


Approved  for  public  release;  distribution  is  unlimited. 

24 
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V0=932m/s 
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Fig.  B-2  Comparison  of  different  values  of  V0 


Changing  projectile  mass  has  very  little  effect  on  oECC  (at  most  between  1%  and 
2%)  (Fig.  B-3).  The  observed  differences  in  aECC  come  from  the  effect  of  changing 
mass  on  lag  time,  which  is  a  primary  driver  of  deflection  due  to  crosswind. 


Fig.  B-3  Comparison  of  different  projectile  masses 

The  subsonic  drag  coefficient  plays  no  role  in  the  trajectory  until  the  range  at  which 
the  projectile  passes  into  the  transonic  drag  regime,  so  Fig.  B-4  shows  a  perfect 
overlap  of  all  3  cases  until  almost  600  m.  The  difference  in  aECC  observed  thereafter 
between  different  cases  does  not  exceed  0.5%. 
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Fig.  B-4  Comparison  of  different  subsonic  drag  coefficients 
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Appendix  C.  Independent  and  Autoregressive  Models 
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C.l  Independent  Measurements 


The  simplest  choice  was  to  draw  the  wind  at  each  station  independently  using  /rn 
and  an.  This  guaranteed  that  the  behavior  of  the  wind  at  each  station  was  as 
observed. 


To  assess  how  well  this  independent  model  performed,  we  calculated  the 
Equivalent  Constant  Crosswind  (ECC)  values  across  the  range  in  question  for  each 
variable  wind  profde  generated  (ECCi),  aggregated  them,  and  computed  their  mean 
and  standard  deviation  to  compare  with  those  of  the  observed  data.  Figure  C-l 
shows  [J.eccb  ar>d  Hecc/  versus  range,  and  Fig.  C-2  shows  <jEccb  and  oECCl  versus 
range. 


Fig.  C-l  Comparison  of  hECCb  and  nECCl 
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Fig.  C-2  Comparison  of  oECCb  and  oECC{ 


Eccb  and  H-ecc,  were  within  0.3%  of  each  other  at  all  ranges,  but  oECCb  and  oECc, 
diverged  with  range,  reaching  a  discrepancy  of  more  than  40%  by  the  maximum 
range.  There  is  much  more  variation  in  ECCb  than  in  ECCi.  This  is  at  least  in  part 
due  to  ECCi  not  taking  into  account  the  correlation  values  between  stations  that  we 
calculated  in  Section  2  of  the  main  report.  Including  those  correlations  called  for  a 
more  sophisticated  model. 


C.2  First-Order  Autoregressive  Model 


An  autoregressive  model  takes  into  account  the  correlation  between  measurements. 
The  basic  fonn  of  a  first-order  autoregressive  model  (so-called  because  it  takes  into 
account  the  first  previous  value  to  generate  the  next  value)  generates  a  sequence  of 
values,  Xn,  each  with  zero  mean  and  unit  variance,  which  depend  upon  the  previous 
states  of  the  sequence,  a  set  of  correlation  coefficients  (pnjrn),  and  a  normally 
distributed  white-noise  process  with  zero  mean  and  unit  variance  (Gn). 1 


Xn  ~  Pn- l.n-^n-l  +  ^7iV^~ 


Pn- 


1,71 


(C-l) 


X1  is  simply  drawn  against  a  normal  distribution  with  zero  mean  and  unit  variance. 
If  we  want  the  members  of  the  sequence  Xn  to  have  a  different  mean  and  standard 
deviation  (denoted  fix  and  ax,  respectively),  we  simply  transfonn  the  sequence  in 
the  following  way: 


'Deserno  M.  How  to  generate  exponentially  correlated  Gaussian  random  numbers.  Los  Angeles  (CA): 
UCLA  Department  of  Chemistry  and  Biochemistry;  2002  Aug  [accessed  2016  Jan  28],  http//www.cmu 
.edu/biolphys/desemo/pdf/corr  _gaussian_random.pdf. 
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(C-2) 


Wn  =  Px  +  axXn 


In  our  case,  the  members  of  the  sequence  represent  the  anemometers  in  ascending 
order.  Since  each  member  of  the  sequence  MVn  has  the  same  mean  and  standard 
deviation  and  our  anemometer  measurements  have  different  means  and  standard 
deviations,  we  need  a  modified  formula  to  preserve  those  measurements.1 

Wn  =  Pn-l,n  (Wn-1  ~  Pn-l)  +  <rnGnyjl  _  Pn-l,n2  +  Pn  ■  (C-3) 

°n- 1 

Deserno1  did  not  include  the  calculations,  but  we  can  easily  show  that  in  this 
formulation,  each  member  of  the  sequence,  Wn,  has  mean  pn  and  standard  deviation 
ern  (variance  ern2).  First  in  the  sequence  is  the  mean. 

E (Wn)  =  E(pn_Xn  (Wn-i  ~  pn-i)  +  ^W1  _  Pn-l,n  +  Pn )  • 

°n- 1 

E(Wn)  =  E (Pn_1  n  (VFn-i  -  Pn- 1))  +  E(anGnJl  -  pn_^)  +  E (pn)  . 

°n- 1 

E (Wn)  =  Pn-l,n-r2^E(Wn_1  -  pn-i_)  +  (Jnj l  ~  pn-i,n2 E (Gn)  +  pn  . 

°n- 1 

E(Wn)  =  pn-1,n  c E(Wn_1 )  -  E(pn_ ,))  +  anJl  -  Pn-l,n2 ' 0  +  Pn  ■ 

°n- 1 

E  (Wn)  —  Pn-l,n  ~  (.Pn— 1  Pn— l)  "F  Pn  ■ 

an- 1 

E(Wn)  =  pn  . 


Next  is  the  variance  (using  the  fact  that  Gn  and  Wn  are  independent,  and  therefore 
have  zero  covariance). 

Var(yvn)  =  F[(pn-i  n  -r2-  (Wn- 1  ~  Pn- 1)  +  (TnGn^l  ~  Pn-l.n  +  Pn  ~  Pn)2]  • 

On- 1 

Var(Wn)  =  E[(pn_ljn2  -^E(Wn^  -  pn_P2  +  an2G2{  1  -  pn-Ul2)  +  2pn_ln^Ljl  -  pn_1,n2(lVn_1  -  pn-PGn\ 

°n- l  °n-l 

Var(Wn)  =  pn_ln2  -  pn.p2  +  (Tn2(l  -  plx-l,n2)E(Gn2)  +  T^V1  ~  />n-U2£((^-l  -  V-n- l)Gn) 


Var(Wn)  =  Pn-I,n2°n2  +  ® n2( 1  “  Pn-l,n2)  +  ^Pn-l.n  f2- a/1  “  Pn-l,n2  ‘  0  . 

°n-l 

Far(fFn)  =  pn_i,n2ern2  +  an2(l  -  p„_ i,n2)  • 

Far(fFn)  =  <7n2  . 


Using  pn,  ern,  and  pn  m  from  Section  2  of  the  report,  we  can  now  implement  the 
first-order  autoregressive  model  and  compare  its  ECC  results  (ECCari)  with  the 
observed  data  (for  all  subsequent  models  we  will  only  show  the  charts  illustrating 
°ecc  because  pECC  was  accurate  to  within  0.5%  at  any  range  using  any  proposed 
model).  This  is  shown  in  Fig.  C-3  along  with  the  results  from  the  independent 
model. 
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Fig.  C-3  Comparison  of  aECCg,  ffECCf,  and  <tECCari 


As  shown,  <tECcar1  diverges  far  less  from  oECCb  across  all  ranges  than  oECCl  but  it 
still  underestimates  <jECCb  by  almost  20%  at  the  longer  ranges.  This  raises  the 
question  of  whether  it  is  valuable  to  take  more  correlations  into  account. 


C.3  Higher-Order  Autoregressive  Models 

Deserno1  only  provides  the  technique  to  implement  a  first-order  autoregressive 
model.  We  derived  a  general  fonnula  to  model  a  given  order  of  autoregression  and 
also  preserve  /rn  and  an.  As  an  example,  we  look  at  extending  from  a  first-order  to 
a  second-order  autoregressive  model.  That  would  involve  an  expression  like  the 
following2: 


vvn  —  Hn- 


But  that  expression,  though  it  takes  into  account  the  previous  2  values  of  Wn  and 
the  corresponding  correlation  coefficients  and  preserves  jun  (this  can  be  verified 
analogously  to  the  process  shown  in  Section  5.2  and  is  not  shown  here),  does  not 
preserve  the  variance  as  ern.2  Finding  the  variance  of  Eq.  C-4  begins  as  it  did  in  the 
first-order  case,  and  the  steps  relying  on  independence  of  the  Gn  terms  with  the 
other  variables  will  be  omitted. 


Var(Wn)  =  E(pn_ln  —s^(Wn_1 

°n-l 


'  pn-l)  "T  Pn-l,n2  Pn-2,n  (Wn-2  P71-2)  ^  ^"rTnr2\T  Pn-2,n2) 


“We  specify  Gn  l  and  Gn  2 


to  clarify  that  the  white  noise  processes  are  separate  and  independent. 


Approved  for  public  release;  distribution  is  unlimited. 

31 


Var(Wn)  ^  (fin-1, n 


2  -^4(W„_i  -  /in_!)2  +  an2GnA2(l  -  pn_lin2)  +pn-2,n  -^4  (Wn-2  ~  Pn-2)2 
°n- l  2  °n-2 

n  Gn,2  (l  ~  Pn-2,n  )  2  Pn-l,nPn-2,n~  “  O^n-1  —  Mn-l)0^n-2  ~  Pn- 2))  • 

°n-l°n-2 

Kar(Wn)  =  2ern2  +  2  ■  pn-i,nPn-2,n  — E((Wn_  x  -  pn-i)  Wi-2  ~  Pn-z))- 


Kar(Wn)  =  2ern2  +  2  ■  pn_ljnpn_ 2, 


,n 


"Pn-Z.n-l^n-l^n-Z 


°n-l°n-2 

Var(Wn )  =  2(7n2  +  2  ■  Pn-i,nPn-2,nPn-2,n-l^n2  • 
Kar(Wn)  =  ern2(2  +  2pn_lnpn_2inpn_2;n_1)  . 


This  points  the  way  forward  for  how  to  modify  our  expression  to  preserve  the 
observed  means  and  standard  deviations  at  the  stations,  which  is  accomplished  in 
the  following  expression: 

Pn-i,n'Z~^(Wn-i~tln-i)+fTnGn,iJ^~Pn-i,n2+Pn-2,n'Z~^(Wn-2~tln-2')+fTnGn,2.P~Pn-2,n2 

wn  =  — - * - +  Mn  .  (C-5) 

More  generally,  we  can  derive  an  expression  that  will  preserve  pn  and  an  for  any 
order  p  of  autoregression.  We  assume  that  for  an  autoregressive  model  of  order  p, 
if  n  <  p  +  1,  the  model  will  use  all  available  previous  measurements  to  generate 
the  nth  measurement. 


Wn=pn+YJ\=1 


( Pn-i,na  .0 ^n-i  Pn-O^^n^n.i  jl  Pn-i.n 2) 
°n-i  \ 


where  L  and  C  are  constants  determined  by 

L  =  min (p,n  —  1)  . 


c  -  A  L  +  U=i  Lj=i+ 1  2 p  n-i,nPn-j,nPn-j,n-i  • 


(C-6) 


(C-7) 

(C-8) 


With  this  expression  in  hand,  we  can  implement  higher-order  autoregressive 
models  and  compare  their  ECC  results  with  the  observed  data.  Fig.  C-4  shows  the 
results  from  using  autoregressive  models  up  to  order  5.  All  curves  represent  ECC 
detennined  using  the  M855/M4  ammunition/weapon  configuration. 
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Fig.  C-4  Multiple-order  autoregression  oECC  results 

Each  higher  order  of  autoregression  produces  results  that  are  marginally  closer  to 
the  observed  data  than  the  next  lower  level,  but  there  is  little  difference  between  the 
models’  perfonnance  past  the  third-order  model.  The  fourth-  and  fifth-order  models 
overlap  almost  completely  and  produce  only  minimal  gains  over  the  third-order 
model.  <JeccAR3  stays  within  10%  of  oECCb  across  the  whole  range,  though  it 
persistently  underestimates  <tECCb  (as  do  all  of  the  models’  results). 

Figure  C-4  and  the  persistent  gap  between  the  models’  results  and  the  observed 
results  raise  the  question  of  whether  the  models  tend  to  accumulate  errors  with 
increasing  range  or  whether  some  other  phenomenon  could  account  for  the  gap. 
One  way  we  can  address  that  question  is  to  reverse  the  order  of  the  anemometers 
and  detennine  <jECCb  if  we  assumed  that  we  were  standing  at  station  Aio  and  looking 
at  station  Aio.  Fig.  C-5  shows  these  results,  which  once  again  assume  the  M855/M4 
configuration. 


Approved  for  public  release;  distribution  is  unlimited. 


33 


Range  (m) 


Fig.  C-5  Reversed  direction  multiple-order  autoregression  oECC  results 

°eccAR4  stays  consistently  within  0.4%  of  ffEccB  across  all  ranges  when  the 
anemometers  are  reversed. 

The  model’s  close  adherence  to  the  observed  data  in  this  case  should  not  be 
overemphasized,  however.  There  still  remains  a  gap  between  the  modeled  results 
and  the  observed  results  in  some  cases,  and  upon  examination,  while  they  often  do 
a  good  job  of  replicating  the  behavior  of  the  wind  data  from  Fort  Benning,  the 
autoregressive  models  of  all  orders  have  drawbacks  that  significantly  limit  the 
reliability  of  their  performance  in  general. 

Because  the  autoregressive  models  implicitly  order  the  wind  measurements  and 
only  incorporate  crosswinds  from  one  direction  when  generating  new 
measurements,  they  are  not  equipped  to  model  changes  in  correlation  that  are 
otherwise  unpredictable  from  the  foregoing  data.  In  Table  2  of  the  main  report,  the 
average  correlation  between  adjacent  stations  from  Ai  to  A5  is  0.25  while  the 
average  correlation  between  adjacent  stations  from  A6  to  A10  is  0.48.  Nothing  about 
the  correlations  between  Ai  and  A5  indicates  such  a  jump  in  correlation  for  the 
subsequent  stations,  and,  looking  in  the  other  direction,  nothing  about  the 
correlations  between  A6  and  A10  indicates  such  a  drop  in  correlation  for  subsequent 
stations.  Since  the  correlations  are  not  inherently  directional,  they  cannot  be  a 
function  of  the  preceding  correlations  (for  whichever  ordering  is  arbitrarily  chosen) 
and  cannot  be  preserved  without  a  model  that  takes  them  all  into  account  from  the 
start. 
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